Zeno effect and ergodicity in finite-time quantum measurements 
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We demonstrate that an attempt to measure a non-local in time quantity, such as the time 
average {A) T of a dynamical variable A, by separating Feynman paths into ever narrower exclusive 
classes traps the system in eigensubspaces of the corresponding operator A. Conversely, in a long 
measurement of {A) T to a finite accuracy, the system explores its Hilbert space and is driven to 
a universal steady state in which von Neumann ensemble average of A coincides with (A) T . Both 
effects are conveniently analysed in terms of singularities and critical points of the corresponding 
amplitude distribution and the Zeno-like behaviour is shown to be a consequence of conservation of 
probability. 

PACS numbers: PACS numbers: 03.65.Ta,Xp,Yz 



Quantum Zcno effect (see, for instance, [T|- [3J and Refs. 
therein) is often associated with the perturbation fre- 
quent projective von Neumann measurements [4] made 
on the observed system. Suppose, for example, that one 
wishes to determine the duration tq a quantum system 
spends in a particular subspace SI of its Hilbert space. 
Checking ever more frequently whether the system is in- 
deed inside SI one eventually destroys the transitions be- 
tween SI and the rest of the Hilbert space. Thus, because 
of the Zeno effect, a continuoulsy observed system pre- 
pared inside SI would spend there all available time, while 
for a system initially outside SI, tq would be exactly zero 
[5]. Alternatively, one can perform a finite time mea- 
surement of tq in which a meter monitors the system 
over a finite period of time, and one single observation 
is made at the end of the run. One example of such a 
meter is Larmor clock consisting of a spin which rotates 
only when the system resides in S7. From the clock's final 
orientation one is able to determine the value of r, but 
can learn nothing about the precise moments the system 
enters and leaves SI. A conceptually similar measure- 
ment of the duration r q a qubit spends in its state \q) 
proposed in [7] employs a large number of bosons (e.g, 
a weakly interacting Bose-Einstein condensate) trapped 
in one of the wells of a symmetric double well potential. 
With atomic current between the wells increased when- 
ever the qubit occupies the state of interest, the number 
of bosons found in the other well contains, like a reading 
of a Larmor clock, information about r(q). The purpose 
of this paper is twofold: we analyse the Zeno effect arising 
in the high accuracy limit of a general finite-time mea- 
surement. We also study the effects of a measurement of 
long-time average of a quantum variable and search for 
any evidence of ergodic behaviour. Various approaches 
to quantum ergodicity can be found in Refs. [5] - 
with the importance of measurement (s) performed on the 
system emphasised in |10) . A variant of the Zeno effect 
arising solely from strong interaction between a system 
and its environment has been studied in [3J. 

Consider a quantum system in an iV-dimensional 



Hilbert space with a Hamiltonian H. Choosing an or- 
thogonal basis \a m ), m = 1,2, ..N in which H is not 
diagonal, we can write a transition amplitude between 
initial and final states \i) and |/) over a time T as a sum 
over Feynman paths 

N 

U^(T) ee (f\ exp(-iHT\i) = lirriK^oo £ (f\a mK ) x 

mi,..,mK 

K-l 

Y\ ((hn j+1 \ exp(-i^e|a TOj )(a TOl |i) = ^ Uf*~ % \path], (1) 

j — l paths 

where e ee T/(K — 1) and each path is defined by a se- 
quence {mi, tyik} numbering the states \a m ) through 
which the system passes until reaching the final state |/). 
Consider a quantity A represented by an operator A di- 
agonal in the chosen representation, 
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where P n are mutually orthogonal projectors on (one- 
or multidimensional) subspaces spanned by vectors \a m ) 
corresponding to the same eigenvalue A n . We wish to 
measure the value of a Feynman functional 
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where /3(t) is a known function [T^] and A(t) denotes the 
(highly irregular) function traced by the value of A along 
a given Feynman path. Equation |3]) may, for example, 
represent the time average of a quantity A if one chooses 
(3(t) = l/T = const and, in particular, the fraction of the 
time, tq/T, the system has spent in if A is also chosen 
to be the projector onto a subspace SI spanned by \a n ), 
A = Pn = J2 n en l a ")( a "l IX- Tne probability amplitude 
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for the value of F to be y is given by the restricted path 
sum 

&^ l (y,T) = ^ 5{F[path]-y)U f ^ l [path]. (4) 

paths 

where 8(z) is the Dirac delta. Without loss of generality 
we choose the measuring device to be a Neumann pointer 
with position y which interacts with the system over a 
time T, the full Hamiltonian being H(t) = H — id y j3(t)A. 
At t = T the pointer states of the meter, \y), are 
entangled with the system's states obtained by prop- 
agation along Feynman paths satisfying the condition 
F[path] = y |13j . In particular, for the system and the 
meter prepared at t = in an product state |i)|G), the 
probability amplitude to find at t = T the system in 
the state |/) and, simultaneously, the pointer reading y, 
\E , ^ <_l (y), is given by 

¥^(y,T) ee (yK/lexpH f %{t)dt]\i)\G) = 
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J G(y-y')^{y',T)dy, (5) 

where G{y) = (y\G). Note that the second of Eqs.([5| 
is not specific to our choice of von Neumann meter and 
occurs for a wider class of quantum measurements, e.g., 
for the one considered in [7|. With G(y) narrowly peaked 
around the origin, Feynman paths with different values 
of F[path] contribute to different final meter states, and 
l*^^ l (j/)| 2 yields the probability to reach the final state 
and obtain the value y for the quantity in Eq.([3]). The 
measured result contains, however, an intrinsic quantum 
uncertainty as the values of y' within the peak's width 
around y remain indistinguishable. With many differ- 
ent Feynman paths contributing to the transition ([I]) one 
might expect the values of F to have a broad distribu- 
tion, which a more accurate measurement would resolve 
in ever greater detail. The accuracy can be improved by 
making the initial pointer state narrower in the coordi- 
nate space, e.g., by replacing G[y) with 

G a (y) = aV*G(ay), Hm Q ^ 00 |G Q (y)| 2 = S(y) (6) 

where the factor a 1 ' 2 ensures the correct normalisation 
of the new state, / \G a {y)\ 2 dy = f \G(y)\ 2 dy = 1 . 

With the help of Eqs.([5]) and ^ we can now prove 
a general result: an accurate finite time measurement, 
a — > oo, T < oo, would indicate that at all times A(t) 
maintains a constant value equal to one of the eigenval- 
ues An, with the value of the functional ^ equal ex- 
actly to A^ Jq fi(t)dt. The proof follows from observing 
first that cannot be a smooth function for all 

final states |/). Indeed, increasing a while maintaining 
unit normalisation will cause the integral J G a (y)dy — 
cT x l 2 J G(y)dy to vanish. With it would also van- 

ish lim™ = $*-*(y)a-V 2 / G(y - 

y')dy' — 0, thus contradicting conservation of probabil- 
ity for the pointer. Thus, must have a singular 



part, which we evaluate by rewriting Eq.Q as a Fourier 
integral 

&^{y) = (7) 
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Further, writing exp{-i J^[H + X/3(t)A]dt} = 
Kmjc-,00 nJLi exp{-*[H + Xf3(jT/K)A}T/K} we 
note that limx-^oo exp{-i[^" + Xj3{jT/K)A]T/K} = 
E n exp{-iXp(jT/K)A n } exp(-iP n HPnT/K)Pn + 
0(A _1 ). Defining a Zeno hamiltonian as (c.f. [3]) 

H z =J2 P ™HP n , (8) 
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we obtain 
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where the last terms vanishes as A — > ±oo, 
lim\^ f00 (f\u{X)\i) — 0(A _1 ). Inserting Eq.(|9| into 
Eq.Q yields 

*'^(y) = <f\ exp{-iH z T)P n \i) x (10) 
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where ^i^i„ t u(y) is the smooth Fourier transform of the 
last term in Eq.Q. With the contributions from the 
smooth term vanishing in the limit a — > oo |14j the 
monitored system is seen to undergo unitary evolution 
with a reduced Hamiltonian P n HP n in the subspaces 
corresponding to each of the distinct eigenvalues A n , 
n = 1,2... As in the case of the Zeno effect caused by 
frequent observations, [T]- [3], an accurate finite-time mea- 
surement suppresses transitions between different sub- 
spaces. Taking trace over the system's variables we find 
the probability distribution for the functional ^ 

Um^W^T) ee \&^{y,T)\ 2 w 

v;%-a„ f p(t)dt)(i\P n \i). (ii) 

In particular, a strongly observed system starting in a 
a state corresponding to a non-degenerate eigenvalue, 
\i) = \o.m), Pm — \o-m)(o-m\, would follow a constant 
Feynman path A(t) = Am, thus having the time aver- 
age of A exactly equal to Am, and spending all available 
time in its initial state, tm = T. This failure to find 
real evidence of the irregular virtual motion suggested in 
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Eq.([T]) by separating Feynman paths into ever narrower 
exclusive classes according to the value of a functional 
([3| constitutes the finite-time Zeno effect and is the first 
result of this paper. 

Next we show that the evidence of the virtual motion 
is recovered if one performes a long measurement of an 
arbitrarily high but fixed accuracy, a < oo and T — > oo. 
For simplicity we consider the time average of A, thus 
choosing (3(t) = l/T — const. Changing variables in 
Eq.0, z = X/T, and using spectral representation for 
the operator H + zA yields 

<P f ^ l (y) = (2n)- 1 T x (12) 
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where (H + zA)\ip n (z)) — £ n (z)\ip n (z)) . The long time 
behaviour of &*~ % {y) is now determined by the critical 
points Zn(y) of the exponent in Eq.(12|, 
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d z £ n (z)z=z* = V- 



(13) 

Evaluating the integrals in Eq.([l2|) by the stationary 
phase method yields N rapidly oscillating contributions 
containing factors exp[iS n (y)T] with the phases given by 
the Legendre transforms of £ n (z), S„(y) = z s n y — £ n (z„). 
The critical points y^ of S n (y) are determined by the 
condition z^(y^) — so that from Eq.(13) we have 
y s n = d z £ n (z)\ z= Q. Calculating the derivatives with the 
help of the perturbation theory and evaluating the inte- 
grals in Eq.([5]) by the stationary phase method, we find 

Zimr->oo* /,i_i (y) = 
^Ga(y-<0„|A|^ n »</|^ n )exp(-iE„r)(0 n |i),(14) 
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where H\<f> n ) = E n \<f) n ). Extension to mixed states is 
straightforward. For an initial state pi, all E n non- 
degenerate and all ((f> n \A\ip n ) distinct, the probability 
distribution of the meter's readings is given by 



Um a ^ 0C1 limT^ c 1 W l (y,T) = 

^2 S (y~ {<Pn,\A\^n)){K\pi\(t>n) , 



(15) 



where the order in which the limits are taken is essen- 
tial. We note that a strongly observed system prepared 
with a known energy Ek, (4>n\Pi\4>n) = ^fcn, explores its 
Hilbert space is such a way that the long-time time aver- 
age of a dynamical variable A is sharply defined, with the 
value equal to the ensemble average in the projection von 
Neumann measurement of the operator A, {4>k\A\<f>k) . In 
particular, as T — > oo, the fraction of time a system pre- 
pared in a pure state \4>k) spends in a subspace SI spanned 
by the subset of eigenvectors {\a m ),m € $1} tends to 
the measure of the subset pn = J2men \( a m\4>k}\ 2 This 
ergodic-like [15] property of bound quantum motion, to 
our knowledge not yet discussed in literature, is the sec- 
ond result of this paper. Further, as seen from Eq.(14), 
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FIG. 1. (colour online) a) Amplitude distribution 

hiV'T) for l*> = I 1 ), I/) = \+),H = oja z , A = 
and ujT = 100; b) $ f s ^ oth (y, T) for uT = 50 (solid), ujT = 25 
(dashed), and ujT — 10 (dot-dashed). Also shown is G a (l— y) 
[X7S) for a = 100 (thick solid). 



a system starting in a mixed state pi, undergoes relax- 
ation to the same steady state diagonal in the energy 
representation, p f = J2n\M((f>n\Pi\^n)(^n\, regardless 
of the choice of A. Thus, finite-time average of a quan- 
tity A, (A) T ee J yW i (y,T)dy, tends, as T — » oo, to the 
von Neumann ensemble average Tr{Apf}. As in the case 
of a frequently observed system |TTj the equivalence be- 
tween the time- and ensemble averages is established in 
the state produced at the end of measurement, pf, rather 
that in the initial state pi. 

We illustrate the above with a simple example, where 
one wishes to measure the time average of the y- 
component of a spin 1/2, A = a y , for a two-level system 
with Hamiltonian H = ua x prepared in its eigenstate 
|1), H\l,2) = ±w|l,2). Choosing in Eq.||) /3 = const 
we then convert to dimensionless variables, T — > ujT, 
y -> y/y , y = 0T, f3 -> l/T. As in Ref.Q, G(y) in 
Eq.|5]) is chosen to be a Gaussian, 



G a (y) = (2/7r) 1 ^(a) 1 / 2 exp(-a 2 y 2 ). 



(16) 



Figure la shows $Jmooth(y)i *A ± > = ±1 ± ). for 
ujT = 100, with the stationary region clearly seen around 
Vi = Vi — (ll^yl 1 ) = 0- Two Zeno peaks at y/yo = ±1 
predicted by Eq.(ll) and a single 'ergodic' (for want 
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FIG. 2. (colour online) Probability distribution of the time 
average of the i/-component of the spin, A = a y , for a two-level 
system with the Hamiltonian H = log x vs. the averaging time 
T. The system is prepared in the eigenstate |1), H\l) = w|l), 
and the accuracy of the measurement is Ay/yo = 0.05. The 
ergodic peak and two Zeno peaks are indicated by white and 
black arrows, respectively. 



of a better word) peak predicted by Eq.(15) are shown 
in Fig. 2. Transition between the regimes described by 
Eqs.(ll) and |l5|) occurs as the Zeno peaks, however 



tributions from ^i^oothiy)' whose oscillations become 
more rapid as the time T increases. This is illustrated in 
Fig. lb, for the spin-1/2 system described above, with 
= 1) given by Eq.(|^). For wT = 10 the con- 
tribution from $tmooth(y) (dot-dashed) is negligible, and 
the Zeno peak is formed by the ^-singularity of $ +<_1 



(c.f. Eq. (10)) at y = 1. For ujT = 50 the contribution of 



narrow G a may be, are eventually cancelled by the con- 



the singular term is largely cancelled by the first negative 
oscillation of ^tmooth (u) (solid) that fits under the Gaus- 
sian G a (l — y) (thick solid), thus making \E ,+< " 1 (y = 1) 
negligible. 

In summary, we have considered a general finite-time 
measurement based on separating Feynman paths into 
exclusive classes according to the value of a functional 
such as time average of a dynamical variable A, (A) T . We 
have shown that: (i) A highly accurate measurement of 
a fixed duration traps the measured system in the eigen- 
states (eigensubspaces) of the corresponding operator A. 
(ii) For any quantity A, a prolonged measurement of an 
arbitrary but fixed accuracy destroys coherences in the 
energy representation, thus leaving the system in a steady 
state with (A) T equal to the von Neumann projection av- 
erage of A. In the special case of a system prepared in 
a pure stationary state, (A) T is sharply defined, and the 
proportion of time spent in a given subspace f2 is exactly 
equal to the von Neumann probability to find the sys- 
tem there. Both effects are readily explained in terms of 
singularities and critical points of the corresponding am- 
plitude distribution. Present analysis can be extended 
to the cases of 'self-measurement', such as wavepacket 
tunnelling |16| . where no external meter is employed. 
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